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ABSTRACT 


In the network interdiction problem, an interdictor destroys a set of arcs in a 
capacitated network through which an adversary will maximize flow. The interdictor’s 
primary objective is to use his limited resources to minimize that maximum flow, but other 
objectives may be important. Therefore, we describe algorithms for enumerating near- 
optimal interdiction sets in planar networks so that these sets may be evaluated with 
respect to secondary criteria, e.g., safety of attacking forces, collateral damage, etc. Our 
algorithms are based on enumerating near-shortest paths or cycles in the dual of a planar 
network; they find a single optimal interdiction set in pseudo-polynomial time. We 
inclement one algorithm applicable to s-t planar networks (s and t must lie on the 
perimeter of the network) and solve problems with up to 800 nodes and 1274 arcs. An 
exanple of conputational results on that largest network is: The algorithm enumerates all 
19 solutions that are within 50% of optimal in 0.15 seconds on a 300 mHz Pentium n PC. 
We also propose, but do not inplement, a somewhat less efficient extension of this 
algorithm to solve problems on general planar networks. 
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EXECUTIVE SUMMARY 


This thesis develops a conputer program to aid in effectively allocating limited 
resources to disrupt the niaximum achievable flow of a single commodity in an 
adversary’s planar, capacitated network. The network might be a transportation or 
pipeline network moving, respectively, war material or oil The interdictor might be 
attacking the network with cruise missiles or aerial sorties. 

Researchers have studied the network interdiction problem over many years, 
motivated originally by the desire to decrease an enemy’s ability to wage war through 
reducing its ability to resupply. More recently, the U.S. government’s attempt to curb the 
flow of illegal drugs into the country has driven advances in network interdiction 
techniques. These techniques focus on finding and implementing an optimal attack 
strategy on the network. However, an interdictor may be faced with a very conplicated 
set of constraints, or at least a long list of desired qualities of an interdiction decision. 

For instance, the interdictor may be concerned about the safety of his attacking forces, or 
unnecessary collateral damage, etc. So, instead of trying to provide the single best 
solution with respect to a long list of criteria, we may prefer to provide the interdictor 
with a set of near-optimal solutions (with respect to some simple criterion such as 
maximum flow), and let the interdictor subsequently evaluate each of the candidate 
solutions in regard to the more conplicated, perhaps even non-quantifiable criteria. 

We address this void in the research by proposing two algorithms that enumerate 
near-optimal strategies. These algorithms build upon previous work in the field of 
network interdiction, and are divided into several sub-algorithms that incorporate typical 
“network techniques.” 
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The first interdiction algorithm is applicable only to planar networks where both 
the sink and soince lie on the outer perimeter of the network. Its &st sub-algorithm 
constructs the dual network from the primal network. The dual is then expanded to deal 
with the consumption of varying, integral amounts of interdiction resource. The second 
solves the interdiction problem by finding a shortest path from the expanded dual source 
node to the closest expanded dual sink node. It does this by working backward from all 
the expanded dual smk nodes. In doing this confutation, it also computes the shortest 
path distances from each expanded node to the closest expanded instance of the sink 
node. Using these distances, the third sub-algorithm enumerates near-optimal paths in 
the expanded dual network; these paths correspond to all of the near-optimal solutions of 
the network interdiction problem on the original network although some duplicates are 
created. These duplicates are easily culled. 

We implement this algorithm and solve problems with up to 800 nodes and 1247 
arcs. Representative conputational results on this largest network are: The algorithm 
enumerates all 19 solutions that are within 50% of optimal in .15 seconds on a 300 mHz 
Pentium n PC. 

The second interdiction algorithm we propose, but do not inplement, is an 
extension of the first. It does require that the network be planar, but the source and sink 
nodes need not lie on the outer perimeter of the networL The algorithm should still be 
reasonably efficient: Computational complexity should increase, in practice, by a factor 
no larger than 2m, where m is the minimum number of arcs among source-sink paths. 




L INTRODUCTION 


This thesis develops a computer program to aid in effectively allocating limited 
resources to disrupt the maximum achievable flow of a single commodity in an adversary’s 
planar, capacitated network. 

A. THE NETWORK INTERDICTION PROBLEM 

In the “network interdiction problem,” an “interdictor” tries to minimize the 
maximum possible flow of a single commodity from a source node to a sink node in a 
capacitated network by efficiently using his limited assets to “interdict,” Le., break 
network arcs; the network user subsequently maximizes flow, subject to arc capacities, 
flow balance constraints and the fact that interdicted arcs are unusable, i.e., have zero 
capacity. Researchers have studied the network interdiction problem over many years, 
motivated originally by the desire to decrease an enemy’s ability to wage war through 
reducing its ability to resupply (WoUmer 1964, Preston 1970). More recently, the U. S. 
government’s attenr 5 )t to curb the flow of illegal drugs into the country has driven 
advances in network interdiction techniques (Steinrauf 1991, Wood 1993). 

The network interdictor is concerned with finding and implementing an optimal 
attack strategy on a network. Models in the literature have defined “optimal” to mean 
minimizing the maximum flow subject to the number of arcs destroyed being limited to a 
specified number, or subject to other simple constraints. However, an interdictor may be 
faced with a very complicated set of constraints, or at least a long list of desired qualities 
of an interdiction decision. Instead of trying to provide the single best solution with 
respect to a long list of criteria, we may prefer to provide the interdictor with a set of 
near-optimal solutions (with respect to some simple criterion such as maximum flow), and 
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let the interdictor subsequently evaluate each of the candidate solutions in regard to the 
more complicated criteria. Specifically, enumeration of near-optimal strategies, which we 
also call ‘partial enumeration,” provides another degree of freedom for choosing a good 
attack strategy. For instance, the interdictor may not want to inclement a strategy that 
isolates a specific node by destroying all adjacent arcs. In addition, the interdictor may be 
concerned about the safety of his attacking forces, or unnecessary collateral damage, etc. 
However, these strategies may be optimal in the simple, mathematical sense. Through 
partial enumeration, a strategy may be revealed that does not directly isolate the node, yet 
is within a few percent of being optimal with respect to a single criterion — and the 
modest loss in “optimality” may not be critical to the interdictor’s objectives. We will 
therefore explore techniques that will allow for the enumeration of near-optimal 
interdiction solutions. 

B. BACKGROUND 

Most work on the interdiction of flow networks is based on the max flow - min cut 
theorem, which states that the maximum flow in a capacitated network is equal to the 
capacity of the minimum capacity cutset (Ford and Fulkerson 1962). The network 
interdiction problem is then a problem of finding an attack strategy that minimizes the 
minimum capacity cut in the after-interdiction network. 

Steinrauf (1991) and Wood (1993) both use mathematical programming 
techniques to solve this network interdiction problem. To illustrate their methods, let G = 
(NA) be a directed network with node set N and arc set A, where an arc is an ordered pair 
(ij) with ij € N. Each arc (ij) has an associated capacity uij and an interdiction resource 
expenditure The decision variables are Xy, the amount of flow the network user will 
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transport along the arc, and jfj, a binary variable equal to 1 if the interdictor chooses to 
interdict the arc, and 0 otherwise. The network interdiction model can then be formulated 
as 

min max 

X 

s.t - X. = 0 

^Xij - ^xji = 0 '^ie N-s-t (1) 

j:(iJ)eA 

^xtj - ^Xjt + xis = 0 

0<Xy <Uij{l - 75) V(i, j)e A 

where F = {75:75e{0,l} V(i\/)e A, Note that by setting '}ij = 0ViiJ) e A, 

the inner maximization reduces to the standard maximum flow model Our basic model is 
identical to this except that we assume arcs are undirected and the network must be 
planar. Wood also shows that by taking the dual of the inner maximization and 
linearizing, the above model can be converted into an integer program so that the problem 

' tf-. 

can be solved using standard integer programming techniques. 

The mathematical programming approach described above is very general For 
instance, it can handle directed, undirected planar, and non-planar networks as well as 
multiple sources, sinks and resource types. However, the resulting models can be hard to 
solve (Cormican, 1995) and adding a partial enumeration capability could be difficult. 
Consequently, we base our techniques on simpler “network algorithms.” TTiese are 
algorithms, like shortest path algorithms, that work directly with network data structures 
rather than linear or integer programming algorithms that operate on an algebraic 
description of the problem. 
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WoUmer (1964) develops a network algorithm that finds the n arcs that, when 
removed firom an undirected planar network, minimize the maximum flow. This algorithm 
makes use of a modified topological dual of the pnmal network (for a detailed discussion 
of the topological dual, see Ahuja, et al. 1993, pp. 262-265) where a new arc of zero 
length is added in parallel to each dual arc. Using a dynamic-programming method, the 
algorithm then calculates the shortest path through the modified dual network, traversing 
at most n arcs of zero length. The zero-length arcs correspond to those that should be 
interdicted in the primal network to minimize the maximum flow. Though WoUmer’s 
algorithm solves the problem in polynomial time, it is limited to an s-t planar network, and 
it assumes that the dual is given and that every arc has an interdiction cost of one unit. 

(An “s-t planar network” has its source node s and sink node t on the periphery of the 
region enclosed by the network; s and t are positioned arbitrarily in a “general planar 
network.”) 

Phillips (1992) removes some of WoUmer’s limitations and solves a variant of the 
network interdiction problem for undirected general planar networks with non-equal arc 
interdiction costs. She describes, but does not implement, three pseudopolynomial-time 
network algorithms that utilize dynamic programming. By introducing an error parameter 
E > 0, she shows how a solution within 1006% of optimal can be returned in polynomial 
time. 

C. PROPOSED ALGORITHMS 

We propose two specialized network algorithms that enumerate near-optimal 
network interdiction solutions. The first algorithm applies to undirected s-t planar 
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networks, and the second to general undirected planar networks. We inclement the frst 
algorithm but not the second. 

Our algorithm for undirected s-t planar networks provides a stand-alone tool that 
takes as its input the adversary’s capacitated transshipment network along with arc 
capacities Uij, interdiction costs r,y, and a description of the embedding of that network in 
the plane. It then constructs the s-t dual network using an algorithm motivated by 
Bhaskar and Sahni (1988) for finding a rectangular dual of a planar triangulated graph. 

Our algorithm then solves for the optimal interdiction using an extension of WoUmer 
(1964). Finally, the algorithm enumerates near-optimal paths in an ejqpanded dual network 
using a dynamic-programming technique similar to Byers and Waterman (1984) for 
finding near-optimal paths in a standard shortest path problem. Some duplicate solutions 
are created, but these are easily culled. 

Our proposed methods are not as easy to generalize as the mathematical 
programming models of Steinrauf and Wood. However, we believe that the better speed 
at which such specialized algorithms as ours are likely to run, and the value of the 
enumerated near-optimal solutions, make them worth consideration. 

D. THESIS OUTLINE 

The remainder of this thesis is organized as follows. Chapter n gives essential 
definitions and notation. Chapter HI develops the algorithm for s-t planar networks. 
Chapter FV shows how to extend the methodology to general planar networks but does 
not inplement this extension. Chapter V investigates the performance of the algorithm of 
Chapter IE on randomly generated networks of various sizes, and Chapter VI provides 
conclusions and recommendations for further research. 
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n. DEFINITIONS AND NOTATION 

A. PRIMAL NETWORK 

The network we wish to interdict is an undirected planar network with a source 
node s and sink node t. We wiQ represent this network by G = (N^), and refer to it as the 
“prinial network.” N is the set of nodes, and A is the set of arcs (ij) which are unordered, 
distinct pairs firomN. We assume that the network is embedded in the plane, that all node 
positions are known, and that all arcs are non-crossing straight line segments with positive 
length. Each arc (iJ) has an associated integer capacity Ujj and an integer interdiction 
“cost” r.y, Le., a resource expenditure required to reduce flow capacity to 0 with 
probability 1. The network interdictor has an integer amount of resource R available for 
interdiction of this network. Note that any arc (iJ) can be denoted “uninterdictable” with 

rij = 

The network user is trying to maximize the flow of a single commodity from s to 
the t in the primal network. (The inner maximization of Equation (1) with all = 0.) This 

model is also valid for a capacitated undirected network if we make the following 
transformation to the undirected network: Replace every undirected arc (iJ) by two 
symmetric directed arcs (iJ) and (/,i), each with capacity Uij, where M,y is the capacity of 
the original undirected arc (e.g., Ahuja, et al. 1993, p. 168). 

B. DUAL NETWORK 

An undirected planar network G = (NA), with a given embedding in the plane, has 
associated with it a unique undirected planar network G* = (N*A*)- G* is referred to as 
“the dual network.” Each dual arc uniquely intersects a primal arc, so that 1A*| = 
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|A|. The number of dual nodes |/V*| is equal to the number of “faces” in the primal 
network. A “face” is a region defined by the arcs in a network, the unbounded region 
being the “outer face” (e.g., Lawler 1976, p. 33). For example, the sequence of nodes 1- 
2-6-5-1 in Figure 1 defines a face of G, and the sequence of nodes 1-2-3-4-8-7-6-5-1 
defines the outer face of G. 

The dual network for a general undirected (primal) network is constructed by 
placing a dual node in each face of the primal network, including the outer face, and then 
connecting dual nodes in adjacent faces with dual arcs. For our purposes, we will also 
assign length Ui*j* = UyBud resource ri*j* = r,y to the dual arc if it intersects primal 
arc (iV). An exanple of this construction is shown in Figure 1. We will refer to this 
network as the “general dual network.” Although the primal network’s arcs can all be 
straight line segments, the dual arcs need not be. 



!-Dual arc 


# Primal node 
O Dual node 


Figure 1. General Dual Network. General planar network 
and its dual Primal arc labels are (Uij,rij). 
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When the primal network is s-t planar, we add an artificial arc from 5 to r to create 
an additional face. The dual network is then constructed by placing the dual source node 
in this additional face, the dual sink node t* in the outer face, and by constructing arcs 
as in the general case but with the omission of the arc This construction is shown 

in Figure 2. We will refer to this network as the dual network.” 


f* 



ArlifidBlAic 


Figure 2. An s*-f* Dual Network. An s-t planar network 
and its dual Primal arc labels are (Uij^ij). 

The general dual network is of interest because of the one-to-one correspondence 
between s-t cutsets in the general primal network and simple cycles with odd parity in the 
dual network (e.g., Phillips 1992). A “cycle” is a path that begins and ends with the same 
node, and a “sirrple cycle” is a cycle which has no repeated nodes except the first. If we 
define P as the set of arcs in any simple path (no nodes may repeat) from ^ to r in the 
primal network, and P* as the set of corresponding arcs in the dual network, “parity” is 
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defined as the number of arcs fi-om the set P* that are traversed in a dual cycle. By the 
Jordan curve theorem, a cycle separates nodes 5 and r if and only if it has odd parity 
(Philhps 1992). Therefore, identification of a minimum capacity s-t cutset in G becomes a 
problem of finding a shortest (in terms of the lengths Ui*j*) simple cycle in G* with odd 
parity. 

For s-t planar networks, a single cycle with odd parity reduces to an s*-t* path in 
the s*-t* dual network (e.g., Ahuja, et al. 1993, p. 263). In terms of the 5 *-r* dual 
network, identification of the minimum capacity s-t cutset becomes the problem of finding 
a shortest path in G* from s* to t* with respect to arc lengths «;»/». 

C. EXPANDED DUAL NETWORK 

We now modify the dual network so that we can use a shortest path algorithm to 
solve the network interdiction problem. We call the modified dual network the “expanded 
dual network” and denote it as G^ = (The expanded dual network is an 

abstraction, and need not be created when actually solving the problem.) In addition, the 
modification varies slightly depending on whether the dual network is an s*-t* dual or a 
general dual We will discuss the expanded s*-t* dual network first and then build upon 
that descr^tion for the expanded general dual network. 

1. Expanded s*-t* Dual Network 

Given that an interdictor has R units of resomce available for interdiction, the 
expansion creates R +1 “expanded nodes” for each dual node. Each expanded node is 
identified by an ordered pair where i* is the original dual node label and 
r € {0,...,/?} is the level of resource expenditure. Starting from any expanded node 
there are two options given an incident arc in the s*-t* dual network: 
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(a) Traverse the “ejq)anded arc” whose length is Ui*j*, to the node {j*,r}; 

this is equivalent to leaving the corresponding primal arc intact; or (b) if r + r,*/* < R, 
traverse the expanded arc + r,*/*}), whose length is 0, to the node {j*,r + 

n*/*}; this is equivalent to interdicting the corresponding primal arc (WoUmer 1964). An 
exanple of a single s*-t* dual network and the corresponding ejq)anded s*-t* dual 
network with /? = 1 is shown in Figure 3. The network'interdiction problem in terms of 
the expanded s*-t* dual network reduces to finding the shortest of the shortest (distance) 
paths fi'om {s*,0} to for re {0,...^}. 
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Figur6 3. Ulustr&ting the exp^ded j*-/* duEl network; (a) a snnple s*-t* dual network 
with arc labels (b) the associated expanded s*-t* dual network for i? = 1 and 

arc labels Ui*j*. 
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2. Expanded General Dual Network 

The expanded general dual network is a sinnple extension of the expanded s*-t* 
dual network. The extension creates two new nodes for every expanded s*-t* dual node, 
one for odd parity and one for even. If an arc e P*, then its parity is 1; otherwise 
an arc’s parity is 0. We then define the parity of a path to be the exclusive-OR of the 
parities of the arcs it contains. An expanded general dual node is identified by an ordered 
triple {i*,r,p} wherep is the (binary) parity. 

Starting from any expanded network node {i*,r,p}, there are four general options 
given an incident arc in the dual network: 

(a) If e P*, traverse the expanded arc ({i*,rp],{j*,r,p}), whose length is m,*/* , to 
the node {y*,r,p}; this is equivalent to leaving intact the corresponding primal arc (ij) 
€P; 

(b) g P* and r + <R, traverse the expanded arc ({i*,r,p},{j*yr + ri*j*,p}) 

whose length is 0 to the node {j*,r + r,-*/* ,p}; this is equivalent to interdicting the 
corresponding primal arc (iJ) ^ P; 

(c) If € P*, and letting p ' denote the parity opposite to p, traverse the expanded 
arc ({i*,r,p},{y*,r,p"}), whose length is , to the node 0; this is equivalent to 
leaving intact the corresponding primal arc (iJ) e P; 

(d) Or, if e P* and r + ri*j*<R, traverse the expanded arc 
i{i*,r,p},{j*,r + ri*j*,p 'j) whose length is 0 to the node {j*,r + r,*/. ,p 'J; this is 
equivalent to interdicting the corresponding primal arc (iJ) € P; 
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The network interdiction problem, for a general planar network, is therefore the 
problem of finding the shortest of the shortest (distance) paths in the expanded general 
dual network fi-om {i*,0,0} to {/*,r,l) over all/* € N* and over allre {0,...,/?} (Phillips 
1992). 
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m. INTERDICTION ALGORITHM FOR s-t PLANAR NETWORKS 


This chapter describes our network interdiction algorithm for to s-t planar 
networks. The algorithm consists of a main program that calls three procedures or “sub¬ 
algorithms.” The first sub-algorithm constructs the dual network from the primal 
network. The second solves the interdiction problem by finding a shortest path from 
expanded node to the closest expanded node of the form by working 
backward from nodes for r e {0,...^}. Simultaneously, it also computes the 
shortest path distances from each expanded node {j*,r} to the closest expanded 
instance of t*. Using the values f(j*,r), the third sub-algorithm enumerates near-optimal 
paths in the expanded dual network; these paths correspond to all of the near-optimal 
solutions of the network interdiction problem on the original network; some duplicate 
interdiction solutions may be generated, but these are easily culled. We first describe the 
data structures used by the sub-algorithms and state assunptions regarding input data. 

A. DATA STRUCTURES 

We assume that the primal network is embedded in the plane so that no arcs cross. 
We assume that all node positions are known, and that all arcs are straight line segments 
with positive length. The input file contains the number of nodes, the Cartesian 
coordinates of each node, and for each arc (ij) the associated values 1,7, %, and rij. 

We utilize three main data structures to implement our interdiction algoritlun. The 
first of these is a network data structure used to store the primal network; this is the 
“primal data structure.” This data structure uses an array of linked lists, one list for each 
primal node. The linked list for each node i contains a record for each adjacent node j 
sorted by increasing 6,y; is defined as the angle in radians measured clockwise from the 
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positive y-axis to the arc (ij) with the Cartesian coordinate system centered at node i. 
Using Figure 1 (in Chapter H), two examples are: Qn = 7 t /2 radians, and Q 26 = it radians. 
Each record in the linked list for node i contains a field for the adjacent node number j, the 
associated arc number, 6y, Uij, r y, and a pointer to the next adjacent node or a null pointer 

if there are no more adjacent nodes. For an exantple of a primal data structure, see Figure 
4. The primal data structure is filled with data fi-om an input file. 



s* 


i 

/ 



Primal arc 

— 

Dual arc 
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Figure 4 The primal data structure. The first node is s and the last is t. (a) an s-t 
(primal) network with arc labels (My,ry) and the associated dual; 

(b) the associated primal data structure for nodes s and t 
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The second data stracture exploits the fact that each arc in the primal network has 
associated with it a unique arc in the dual network; this structiure is the “arc data 
structure.” Its purpose is to allow quick output of a solution in terms of the primal 
network despite the fact that the solution was obtained by manipulation of the dual 
network. The arc data structure consists of an array of records with one record for each 
primal arc (zj), each record containing a field for i,j, Uij, nj, and the dual nodes i*,j*, 
which correspond to the crossing dual arc The first four fields are filled with 

primal network data fi'om the input file, and the last two are filled as the dual network is 
created. 

The third data structure is identical to the first except that it is filled by reading the 
dual network from the arc data structure detailed above once the dual has been created by 
our dual-finding algorithm; this structoe is the “dual data structure.” Although we do not 
create an explicit representation of the expanded dual network, we do need to store node 
information for that network, namely, the values/(/*,r). This is done in a two-dimensional 
array with R [A/] elements. 

B. DUAL-FINDING ALGORITHM 

The first sub-algorithm in our network interdiction algorithm constructs the dual 
network. Discussion of the topological dual of a planar network is ubiquitous in the 
network literature (e.g., Ahuja, et al., 1993, pp. 262-265, Lawler, 1976, pp. 32-36). 
However, we are unable to locate in the literature an algorithm that constructs the dual of 
a planar network. Bhasker and Sahni (1988) do describe an algorithm to construct the 
dual of a planar triangulated network. Unfortunately, the triangulated assunqption is too 
restrictive for our work as it requires: Every face is a triangle, all internal nodes have 
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degree > 4, and all cycles that are not faces have length > 4. However, their algorithm 
does provide a starting point, and given this starting point and O’Rourke’s (1994) brief 
discussion, we are able to derive an algorithm to construct the dual of an s-t planar 
network. 

Given an s-t planar network (as in Figure 5(a)), it is easy to construct the dual (Figure 
5(b)) using the techniques discussed in Chapter n. 



Figure 5. An s-t primal netwoik (a) and 
its dual (b) with arc ideittifiers aixi 
arc data/(^y,r^). 


The key element in constructing the dual is to recognize a face in the primal network. As 
each inner face corresponds to a cycle containing no subcycles, the algorithm must identify 
these special cycles. Our algorithm operates as follows. 

Since each arc lies on the perimeter of two faces, a new directed network G'is 
constructed from G by replacing each arc (ij) with two directed arcs in anti-parallel, (ij) 
and (/,0. In addition, we create two artificial directed arcs, is,t) and (r, 5 ). These arcs 
represent the artificial arc (^,0 that wraps around G and is used in creating the s*-t* dual. 
See Figure 2 in Chapter H. Starting at s, a directed path is then traversed in G'under the 
constraint that doubling back is disallowed, ie., the arc (J,i) can not be added to the path 
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immediately after traversing arc (ij). When an arc (ij) is traversed to reach a new node j, 
the next arc traversed (j^) must form the sharpest possible angle, in a clockwise direction, 
with (iJ). (At a junction, always take the sharpest turn to the right that is available.) 

When some node i is repeated in the path, a special cycle, ie., a face, has been identified. 
All arcs in that cycle are deleted, and the search is continued from i. If i has no incident 
arcs remaining, a node with incident arcs is identified and the search is restarted from 
there. Once all arcs in the network have been deleted, all faces have been identified. 

To implement the algorithm, we utilize the primal data structure which stores each arc 
twice. We do not actually wait until a cycle has been identified to delete arcs, but rather, 
we delete them as they are traversed. (The current path is stored in a stack, so nothing is 
lost by immediately deleting arcs.) The algorithm continues to look for faces until no arcs 
remain. For exanple, the dual of the primal in Figure 5(a) would be found as follows. 
Beginning a.ts(s = 1), the algorithm first traverses to node 2, then to node t, then to node 
3 and then back to s. As this completes a path from j to 5, a face is found which encloses 
dual node 2*. The algorithm then looks for the next face. Since 5 has untraversed 
incident arcs, and since it is where the cycle was identified, the search restarts from there,- 
traverses to node 5, then to node t, then to node 4, and back to node s locating the second 
face, i.e. dual node 5*. The algorithm continues to return to node s until all incident arcs 
have been deleted from it. The algorithm then increments to the next node j = 2 until all 
incident arcs have been deleted from it, then to j = 3, etc. When no primal arcs remain, all 
the dual network nodes have been located. 

The dual construction is conpleted by creating the dual arcs in a process that exploits 
the one-to-one correspondence between primal and dual arcs. We utilize the arc data 
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structure which stores ij, uij, nj, i*, andy* for each of these arc pairs. The first four items 
we assume are given, the last two we define as the dual nodes are found. For example, 
when we have identified a face (which corresponds to dual node k*) with primal arcs (zo, 
h), O’l, hX («Vi, j'o), we set i* = k*, orj* = ^* if i* is already ^ecified, for each arc in 
the path identifying the face. 

1. Assumptions 

The algorithm closely follows the discussion above, but there are a few details that 
should be expanded. In order to construct the dual network, it is vital for us to know the 
location of both s and t in G. Without this information we would be unable to construct 
the artificial arc which wraps around G and connects s to r. By definition of an s-t planar 
network, both s and t lie on the outer perimeter of the network. Since s and t Ue on the 
outer face of G, there exist angles 0,, and 0„ respectively, such that a line segment with 

origin at s (at t) and extending at an angle of 0^ (0„) does not intersect G. Clearly, these 
line segments can be joined into an artificial arc (s,t) which bends around G without 
intersecting it. Our algorithm assumes that 0^ and 0«, as described, are given as inputs. 

The pseudo-code for the algorithm follows. It uses ^ijt s Qjn - Qij. (For exan 5 )le, 
^126 is 7 c/ 2 radians in Figure 1 in Chapter H.) 
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2. Pseudo-Code for Dual-Finding Algoritiim 


Procedure: FINDDUAL(G) 

Input: An undirected planar network G = (N,A) in adjacency-list form. 

(Arcs incident to each node are sorted by increasing 6y.) 

Arc[A] /* an array of records containing i,J, Uij, rij */ 

Output: G* = (N*A*), the dual network associated with G 
in adjacency list form; 

{ 

Add artificial arcs (s,f) and (t,s) to G; 

Mark all f e as “oficycle”; 

DualNodeCounter 0; 

For m = 1 to { 
i m; 

Push i on top of stack S; 

Mark i as “oncycle”; 

While i still has an incident arc { 

Select an arc (ij) which maximizes (|) s.t. radians; 

If no such arc exists, choose the first arc in adjacency list; 
Delete arc (ij) fi'om the adjacency list; 

If j “oncycle” /* face found */ { 

DualNodeCounter= DualNodeCounter+1 
Pop all nodes fi'om S until top of S' = j; 

Mark all nodes popped as “offcycle”; 

For each arc (ij) that was “oncycle” { 

Arc[(iJ)].i* 4- DualNodeCounter; 

} 

i i 

} 

else { 

Push j on top of stack 5; 

Mark j as “oncycle”; 

I ^j; 

} 

}/*j has no more adjacent nodes */ 

Pop / from S; 

Mark i “offcycle”; 

} 

} 
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3. Complexity Analysis 

Each arc (ij) is “selected” once, and this selection process, assuming bounded 
degree, takes 0(|1|) time. Therefore, there is 0(|A|) total work associated with selecting 
arcs. Once an arc is selected, there is only 0(|1|) additional work, so the complexity of 
this algorithm is 0(\A[). Without the assumption of bounded degree, the conplexity is 
^>(|iVIlA|) 

C. DUKSTRA’S SHORTEST PATH ALGORITHM 

The second sub-algorithm now computes the optimal solution to the interdiction 
problem. This solution corresponds to the shortest of the paths in the expanded dual 
network from the expanded node {s*yO} to nodes of the form We accomplish this 
by working backward in an implicit representation of the ejpanded dual network using a 
binary-heap implementation of Dijkstra’s shortest path algorithm (e.g., Ahuja, et al. 1993, 
pp. 115-116). The shortest path distance from any instance t* to itself is 0, so the “reverse 
shortest path distance”/(?*,r) is initialized to 0 for aU r, and the heap is initialized with all 
of these expanded nodes as starting nodes. In computing the optimal interdiction solution, 
the algorithm also computes the shortest path length from every expanded dual node 
{/*,/•} to the closest expanded instance of t*. We denote these distance as/(/*,r), and the 
shortest path length (optimal solution value) asf* =f(s*,0). 

Importantly, the shortest path distances/(/*,r) are exactly what is needed for the 
partial-enumeration sub-algorithm. 

The algorithm proceeds in an obvious manner except that the existence of 
expanded arcs is checked on the fly: Suppose that the algorithm has just removed 
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expanded node {/*,/"} from the heap and must scan arcs incident to this node. For each 
arc in the regular dual network, the expanded dual arc with length 

Ui*j* will typically exist and may be scanned since it corresponds to not interdicting some 
primal arc (ij) and not consuming any resource. (An exception occurs when i*= s* and r 
> 0.) On the other hand, the expanded arc ({i*,r -r^j* },{/*,/■}) cannot exist and will not 
be scanned if r - r,*/* < 0, or r - ri*j* > 0 and i* = s*. 

The corrplexity of a binary-heap inplementation of Dijkstra’s algorithm is 
OdAjloglA^) (e.g., Ahuja, et al. (1993), p. 116). Our modification essentially runs on the 
expanded dual network which has R\N\ nodes and 2R\A\ arcs. Therefore the complexity of 
our algorithm is pseudo-polynomiaL We do not include pseudo-code for the algorithm 
since it is straightforward. 

D. PARTIAL-ENUMERATION ALGORITHM 

The third sub-algorithm is based on Byers and Waterman (1984). Their algorithm 
enumerates all s-t paths within a prescribed distance of the optimal path length from 
somce to sink in a directed acyclic network. They modify a (total) path-enumeration 
algorithm (which is just a modified depth-first search) to locate these paths, and a push¬ 
down stack to store all current valid paths. (The algorithm is related to an “A* search,” 
e.g.. Hart, et al., 1968.) We modify this algorithm to work on an implicit representation 
of the e)q)anded dual network. 

The partial-enumeration sub-algorithm finds near-optimal solutions by traversing all 
paths from {^*,0} to {t*,r},for r e {0,...,R} that are within 100e% of optimal Let 
denote any sirtple path from {5*,0} to {i*,r}, define d{S[i*^)) as its length, and define 
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3S the totdl resource expended to traverse R units of resource are assumed 
available for interdiction. 

Given a partial path the algorithm determines if 

d (5(fv}) + Ui*j* +f(j*,r) < (1 + e)f*. (2) 

If it is, the algorithm feasibly extends the path 5(,*,r) by adding the uninterdicted arc 
and continues the search from with partial path S{j*^r)~ This is equivalent to 
traversing the arc in the expanded dual network. If it is not, the algorithm 

determines if 

^(V.r}) + r,*;^</?and, (3) 

diSn*^)) +/(/*,r + ri*f) < (1 + e)/*. (4) 

If those conditions are met, the algorithm extends by adding the interdicted arc and 
continues the search from {J*,r + rpj*} with partial path }• This is equivalent to 

traversing the arc ({i*,r},{/*,r+ r,-.;.}) in the expanded dual network. The algorithm 
continues this recursion until all near-optimal paths have been enumerated. This method 
can produce the same interdiction set more than once because it enumerates s-t cuts in G 
that may not be minimum. Consequently, duplicate, non-minimum solutions must be 
culled. 

As an example, consider the network illustrated in Rgure 5, when we set /? = 3 and 
£ ~ We obtain an optimal solution of value 100, two near-optimal solutions and three 
duplicates; all six “solutions” are listed in Table 1. 
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Arcs Crossing Optimal and 
Near-Optimal Cutsets 

Residual s-t 
Maximum Flow 

(1,4)' (3,5) (2,5)' 

100 

(1,4)' (1,3)' (2,5) 

102 

(1,4)' (3,5)' (2,5) 

102 

(1,4)' (1,3)' (1,2) dupUcate 

103 

(1,4)' (3,5)' (1,2) duplicate 

103 

(1,4)' (1,3) (2,5)'duplicate 

106 


Table 1. Solution for primal network shown in Figure 5 
with s = l, t = 5,l? = 2 and e = .06. The notation (i,/) indicates 
that the arc is not interdicted; and (i,/)' indicates that it is. 

We demonstrate how the last solution was completed assuming that 5{,»,r) = {{^*,0}, 


{5*,7}}, so that d(S{i*,r]) = 0 and/?(S{iv)) = 1: 

(a) The algorithm will traverse to {2*,7} via the uninterdicted arc B if formula (2) is true. 
As 0 + 106 + 0 < (1 + .06)100 the algorithm does add the node {2*,7} to the path. 

(b) The algorithm will traverse to {t*,l } via the uninterdicted arc A if formula (2) is true. 
As 106 + 103 + 0 > (1 + .06)100 the algorithm cannot add the node {t*,7} to the 
path. 

(c) The algorithm will traverse to {t*,4} via the interdicted arc A if (3) and (4) are true. 
As 1 + 3 > 3 the algorithm cannot add the node {r*,4} to the path. 

(d) The algorithm will traverse to {r*,7} via the uninterdicted arc 7) if formula (2) is true. 
As 106 + 102 + 0 > (1 + .06)100 the algorithm cannot add the node {t*,l } to the 
path. 

(e) The algorithm will traverse to {r*,5} via the interdicted arc D if both formulas (3) and 
(4) are true. As 1 + 2 < 3 and 106 + 0 + 0 < (1 + .06)100 the algorithm does add the 
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node {?*3} to the path. Because this is a sink node, the algorithm outputs the path in 
terms of the primal network as shown in the last row of Table 1. The path is shown in 
terms of the dual in Figure 6. 

(f) Because every arc incident to [2*,1 } has been investigated, the algorithm now 
backtracks to node {5 *,7} and continues investigating uninvestigated arcs. 



Figure 6. Illustrating a solution 
as a dual path with arc labels L 


Let 5 be a stack that stores pairs <j*J}> where j* is a node on the path being 
enumerated and b indicates the status of the arc, say that reaches j* on the current 
path from s*; b = \ inplies that the corresponding primal arc iij) is interdicted; otherwise 
b-0. The pseudo-code for the algorithm follows. 
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1. Pseudo-Code 


Procedure: ENl]MERATE(G*,/,/'', e) 

Input: Dual network G *= (N*A*) in adjacency list form, V /*e N,J* and e; 

Output: All solutions to the network interdiction problem with capacity (1 + e)p; 

{ 

Mark all i*e N* - s* as, “off path”; Mark s* as “on path”; 

Mark all e A* as “untraversed”; 

S 0; 

Push <s*,-l> on top of stack 5;/* “-1” indicates undefined */ 

/* (S, will be an encoding of the dual path from s* to t* along with the status of the 

arcs on that path */ 

^ 0; 

While59£0 { 


} 


<i*,b> <-topof,S; 

While not all arcs out of /* are “traversed” { 

Select an “untraversed” arc 
If;* is “off path” { 

If d (S{iv)) + “<7* + - (1 + £)/*' { 

/* arc is not interdicted */; 

If ;•* = t* { Print(S, }; 

else { Push <j*,0> on 5; 

Mark j* as “on path”; 
i* 

} 

} 

else { 

If ^ R [ 

If d (S„vi) + < (1 + e)/* { 

/* arc is interdicted */ 

If ;•* = r* {Print(S, <t*,l >)}; 
else { Push <j*,l> on S; 

Mark j* as “on path”; 
i* <-;*; 

} 

} 

} 

} 

If arc has been interdicted { Mark (i*,;*) as “traversed”; } 

} 

} 

/* aU arcs out of i* have been “traversed” */ 

Mark aU arcs out of i* as “untraversed”; 

Mark i* as “off path”; 

Pop <i*,b> fi’om S; 
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2. Complexity Analysis 

It is easy to construct a planar network so that its dual has 0(2"^') paths between 
two nodes s* and /*. If most of these paths correspond to candidates for near-optimal 
solutions to the network interdiction problem in the primal network (and this can be 
arranged), then any algorithm to enumerate these paths must be exponential. As each path 
in the expanded dual network might take C>(/?L4I) work to find, the complexity of the 
partial-enumeration algorithm could be 0(i?|A|2'^). In practice, the run time is strongly 
dependent on e and the number of paths enumerated. If e is not too large, and only a few 
near-optimal paths exist, the algorithm is quite efficient. In practice e could be 

incrementally increased firom a very small value until the number of near-optimal paths 
becomes suitably large, but still manageable. 

We perform culling of duplicate solutions off-line using a spreadsheet sort and do 
not include it complexity. Clearly, this procedure could be performed efficiently using a 
has function and table as candidate solutions are generated. 

In Chapter V we will provide con^utational results for the interdiction algorithm 
described in this chapter. Before we do that, however, we describe an algorithm to solve 
the interdiction problem on non s-t planar networks. 
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IV. ALGORITHM FOR GENERAL PLANAR NETWORKS 


This chapter describes, but does not implement, our network interdiction algorithm 
applicable to general s-t planar networks (Le., planar networks where s and t do not 
necessarily lie on the perimeter of the network). We propose foiu sub-algorithms: The 
first conducts a breadth-first search to identify a simple path P between s and t in the 
primal network. The second constructs the (general) dual network from the primal 
network. The third sub-algorithm, a variant of Dijkstra’s algorithm, operates on an 
expanded dual that incorporates parity, Le., the number of times that a path crosses P. To 
solve the interdiction problem, this sub-algorithm must be called once for each dual arc 

in a special subset of dual arcs. A set of restricted optimal \alues is 
computed in this way so that the global optimum isf* = min(,.^)/*'(,.j»). The fourth sub¬ 
algorithm performs enumeration much as in the previous chapter, but must be called for 
each special arc with/*(,. <(l+e)f*. Also, before each enumeration call, the third 
sub-algorithm must be rerun to compute/(,.j») that denotes the length of the 

shortest path firom expanded node {j*,r,p} (p is parity) to the closest node {j*,r,l} for 
re{0,...,/?}. 

A. BREADTH-FIRST SEARCH 

The identification of a minimum capacity s-t cutset in a non-s-r-planar network is 
equivalent to finding the shortest simple cycle in the dual network with odd parity. Here, 
parity counts the number of times the cycle crosses a specified, simple s-t path. (See 
Chapter H.) Our interdiction algorithm is similarly dependent on finding cycles of odd 
parity. Any simple path between s and t will do for counting parity, but for efficiency’s 
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sake we prefer a path with a minimum number of arcs. We identify such a path in 0(|A|) 
time using a standard breadth-first search (e.g., Cormen etal., 1992, pp. 469-475), 

B. DUAL-FINDING ALGORITHM 

Constructing the dual of a general planar network is identical to constructing the 
s*-l* dual as described in Chapter m except: No artificial arc {s,t) is added to G and thus, 
the directed arcs {s,() and (r,^) are not added to G' (the directed network used in 
constructing the dual G*). 

C. A MODHFIED DIJKSTRA’S SHORTEST PATH-ALGORITHM 

The third sub-algorithm, a variant of Dijkstra’s shortest path method, will be cahed 
multiple times to solve the interdiction problem and to help enumerate near-optimal 
solutions. The optimal solution to the interdiction problem corresponds to the shortest 
possible path firom some expanded node {i*,0,0} to a node of the form {i*,r,l} for 
re {0,.. We could call the shortest path algorithm for each i* e N* to solve the 
interdiction problem optimally, but the search can be limited further, as described next. 

An optimal interdiction solution corresponds to a cycle in the dual that crosses 
path P an odd number of times. We may view the optimal cycle as starting at some dual 
node I*, immediately traversing dual arc that crosses arc and then 

eventually returning to i*. TTius, we could solve the interdiction problem by examining all 
shortest cycles that immediately cross from one side of P to the other: 


30 



For each dual arc crossing P: 

Find a shortest path from {i*,0,0} that traverses an expanded 
instance of (i*,j*) to a node of the form } and eventually 
returns to a node of the form {i*,r,l} for re {0,...^}. 

The global optimum corresponds to the shortest shortest path found. Let A** 
denote the special set of arcs from which we initiate the modified shortest path algorithm, 
as above. We would normally expect that L4**I«IA^*I, so significant efficiency should be 
gained by limiting the shortest path searches as described. 

The shortest path algorithm must also be modified to consider only certain arcs 
leaving the source node which is {i*,0,0}in this case. This modification is trivial: 
From{i*,0,0}, scan only “allowable” arcs whose head nodes are of the form }• 
Actually, as in the s-t planar case, the shortest path algorithm will be working backwards 
in the expanded network to reach {i*,0,0}. Therefore, the shortest path search is started 
simultaneously from nodes of the form {i*,r,l} and initially scans arcs entering from nodes 
only. 

Our choice for a shortest path algorithm is the modified version of Dijkstra’s 
method that operates on the (implicit) expanded dual network, as in Chapter HI, but is 
modified further. The algorithm must be adjusted to compute parity properly and to 
implement the special scanning rules described above. The scanning rules are 
straightforward to irrplement; handling parity needs more explanation. 

Like our original shortest path algorithm, this one works backwards, in this case 
from nodes of the form {f*,r,l}for re {0,...,/?} to node{i*,0,0}. The algorithm minimizes 
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/(»v) ,p) which is defined as the shortest path distance fi'om expanded node {h*,r,p} 

to the closest node of the form ). (We use h* and k* as a generic nodes here to 
avoid confusion with the nodes i* and j* used in the previous paragraph.) The algorithm 
handles parity as follows: Suppose the algorithm is scanning the two expanded arcs, which 
correspond to a dual arc (h*,k*), entering the expanded node {k*,r,p}. If does 

not cross the path P, the algorithm determines if (.h*,r,p) and/(,.^) {h*,r-ffi*j^,p) are 

to be updated. (We have scanned the inplicit expanded arcs {{h*,rp},{k*^,p}) and ({/z*, 
r-r^^*J^,p}, If (A*,A*) crosses Z’, the algorithm determines if/(..j., (A* r,l-p) 

(A*,r-r^*j^,l—p) are to be updated. (We have scanned the implicit expanded 
arcs ({A*,r,l-p}, {k*,r,p}) and {{h*, r-rh* 1 ^,\-p}, 

We denote the optimal solution value determined for arc e A**, as above, 

and let/* = min(,.j.)/*(,.j,). So,/* is the optimal solution value to the network 
interdiction problem on G. We will call the shortest path algorithm again, multiple times, 
fi:om within the enumeration algorithm described next. 

D. PARTIAL-ENUMERATION ALGORITHM 

The fourth sub-algorithm differs fi'om the enumeration algorithm detailed in 
Chapter m in that it must be called for each arc (/*/*)€ A**. For each such node, it 
begins by calling the third sub-algorithm to compute the values/;,.^*) {h*,r,p). Once these 
values are obtained, the algorithm searches the expanded network — arcs are generated on 
the fly, of course - looking for paths fi:om {i*,0,0} to {i*,r,/} that are within 100e% of 
optimality. This procedure is identical to the partial-enumeration algorithm in Chapter m 
except that it must account for node parity, just as the shortest path algorithm does. 
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We conjecture that the work performed by both the shortest path and partial- 
enumeration sub-algorithms might be reduced if those procedures could be started from a 
set of arcs that is a proper subset of A**. Suppose we have just identified an optimal 
cycle that starts at (h*,ji*) e A** and then recrosses P twice using ( 12 *,/z*) e A** and 

(.h*,h*) e A**, say. (It must cross P an odd number of times.) Then, we need not begin 
a search for an optimal cycle starting with (h*,]!*) or {h*,jz*) since the cycle we just 
described would just be found again. Thus, the number of calls of Dijkstra’s sub¬ 
algorithm for computing/*' could definitely be reduced. In this solution, it may also be 
possible to reduce the number of calls to the enumeration sub-algorithm, but this is less 
certain and needs further investigation. 

Culling of candidate near-optimal solutions is performed the same way it is for s-t 
planar networks. 
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V. COMPUTATIONAL RESULTS FOR THE INTERDICTION ALGORITHM 


In this chapter we demonstrate the operation of our interdiction algorithm for s-t 
planar networks on four randomly generated problems. The purpose is to demonstrate 
accuracy and efficiency, and potential usefulness to an analyst in choosing a network 
interdiction strategy. We run the algorithm on a small network to illustrate what an 
enumeration output looks like and to allow the reader to verify correctness. We run the 
algorithm on three large networks of varying size to demonstrate computational efficiency, 
and thus the potential for solving large, real-world networks. 

A. ALGORITHM IMPLEMENTATION AND NETWORK GENERATION 

Our algorithm is coded and conpiled using the Borland C-h- Builder 3 and run on 
a 300 mHz Pentium n PC. We also devise a C++ program to randomly generate a planar 
grid network. For an input value «, the program generates an n by 2n grid network. All 
arcs between adjacent nodes on the perimeter of the network are added; this defines the 
outer face. In addition, all other nodes that are adjacent horizontally are connected by 
arcs. The program then completes the network by randomly constructing about 60% of 
the vertical arcs in the interior. The program randomly assigns r.y and % uniformly in 
ranges specified by the user. It then writes this data to an output file that the interdiction 
algorithm can read. An exanple of a randomly generated network for n = 3, Uij e 
{95,...,104}, and r,y e {1,2} is shown in Figure 7; this is also the network we will refer to 
as the “small network.” 
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(99.2) 


( 102 , 1 )^ 


(95,2) 
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(96,2 j (100,2) 


( 102 , 2 ) 




Figure 7. Small Network. The priinal network 
generated forn = 3, s = 1, r = 18, and arc data (My,ry). 

B. COMPUTATIONAL RESULTS FOR THE SMALL NETWORK 

Utilizing the primal network depicted in Figure 7, we run our interdiction 


algorithm for R — 2 and e — .05, and R — 2 and 8 =. 10, where 5 is node 1, and t is node 


18. The time to run both problems is negligible. The output from the algorithm is shown 
in Tables 2 and 3. (No duplicates are generated in this case.) 


Arcs Crossing Optimal and 
Near-Optimal Cutsets 

Residual 1-18 
Maximum Flow 

(1,7)' (1,2) 

95 

(13,14)' (7,8)' (1,2) 

95 

(7,13)' (7, 8)' (1,2) 

95 

(1,7) (1,2)' 

99 

(15,16) (9,10)' (2,3)' 

99 

(15,16) (9,10)' (3,4)' 

99 


Table 2. Output for the small network with R = 2 and e = .05. 

The notation (ij) indicates that the arc was not interdicted, 
and (ijT indicates that the arc was interdicted. There are six 
interdictions within 10% of the maximum interdiction. 
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Arcs Crossing Optimal and 
Near-Optimal Cutsets 

Residual 1-18 
Maximum Flow 

(1,7/ (1, 2) 

95 

(13,14/ (7,8/ (1,2) 

95 

(7,13)' (7, 8/ (1,2) 

95 

(1,7) (1,2/ 

99 

(15,16) (9,10/ (2,3)' 

99 

(15,16) (9,10/ (3,4/ 

99 

(7,13/ (7,8) (1,2/ 

101 

(13,14/ (7,8) (1,2/ 

101 

(17,18)’ (12,18) 

102 

(13,14) (7,8/ (1,2/ 

103 

(14,15) (9,10)' (2,3)'duplicate 

103 

(14,15) (9,10/ (3,4)'dupUcate 

103 

(7,13) (7,8)' (1,2/ dupUcate 

104 

(17,18) (12,18/ 

104 


Table 3. Output for the small network with R = 2 and e = .10. 

The notation {ij) indicates that the arc was not interdicted, 
and iijy indicates that the arc was interdicted. 

The potential usefulness of our algorithm can quickly be seen in the results shown 
in Table 2. There are three optimal solutions shown. However, had we not enumerated 
near-optimal solutions, only one of these would have been found. In addition, through 
partial-enumeration we identify three near-optimal solutions for e = .05, and 8 near- 
optimal solutions for e = .10. This could provide an interdictor with a wide choice of 
good strategies strategies to enploy that also satisfy secondary or tertiary objectives. 

C. COMPUTATIONAL RESULTS FOR THE LARGE NETWORKS 

Here we generate three large networks and run our algorithm on each. We choose Uij 
€ {40,. ..,60} and rij e {1,2} for each of these networks. We also choose to make all 
vertical arcs, and three quarters of the horizontal arcs on the outer face uninterdictable 
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(nj - co). This allows us to run our algorithm for higher values ofR without completely 
disconnecting s and t. For the first large network (LNl), we set n = 10, which creates a 
10 X 20 network with 200 nodes and 315 arcs. The second large network (LN2) is 

generated by setting n = 15, which creates a 15 x 30 network with 450 nodes and 700 
arcs. For the last large network we set n = 20, which generates a 20 x 40 grid network 
with 800 nodes and 1247 arcs. We run the algorithm for R = 3,5 and 10 for all problems. 

Conputational results are summarized in Tables 4, 5 and 6. We do not include 
results for the dual-finding algorithm or Dijkstra’s algorithm because those run times are 
negligible, less than 0.06 seconds in any case. Off-line culling of duplicate solutions is not 
accounted for in the run time either; there are typically very few duplicates generated. 

The number of near-optimal solutions excludes one that is considered “the optimal 
solution.” 



Enumeration 

(secs) 

Total Time 
(secs) 

Near-optimal 

Solutions 

8 

R=3 

lt=5 

R=10 

R=3 

R=5 

R=10 

R=3 

R=5 

R=:10 

.05 

.00 

.00 

.00 

.02 

.02 

.03 

0 

0 

0 

.10 

.00 

.00 

.05 

.02 

.02 

.08 

0 

0 

18 

.20 

.00 

.00 

.28 

.02 

.02 

.31 

0 

0 

135 

.50 

.00 

.00 

.44 

.03 

.02 

.47 

1 

1 

227 

1.00 

.01 

.01 

.61 

.03 

.04 

.64 

4 

9 

323 


Table 4. Computational Results for LNl. Solving the network interdiction 
problem within 20% of optimal on network LNl, with 10 units of interdiction 
resource, enumerates 135 near-optimal solutions in 0.31 seconds. Duplicate 
solutions have been eliminated. 
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Enumeration 

(secs) 

Total Time 
(secs) 

Near-optimal 

Solutions 

8 

R=3 

R=5 

R=10 

R=3 

R:=5 

R=10 

R=3 

R=5 

RrilO 

.05 

.00 

.00 

.35 

.05 

.07 

.42 

0 

5 

239 

.10 

.00 

.00 

.35 

.05 

.07 

.42 

0 

5 

239 

.20 

.00 

.00 

.35 

.05 

.07 

.42 

0 

5 

239 

.50 

.00 

.01 

1.00 

.05 

.07 

1.06 

1 

12 

649 

1.00 

.02 

.04 

1.33 

.07 

.09 

1.39 

9 

39 

855 


Table 5. Computational Results for LN2.. Solving the network interdiction 
problem within 100% of optimal on network LN2, with 5 units of interdiction 
resource, enumerates 39 near-optimal solutions in 0.09 seconds. 



Enumeration 

(secs) 

Total Time 
(secs) 

Near-optimal 

Solutions 

8 

R=3 

R=5 

R=10 

R=3 

R=5 

R=10 

R=3 

R=5 

R=10 

.05 

.01 

.01 

.03 

.10 

.10 

.15 

1 

1 

17 

.10 

.01 

.01 

.03 

.10 

.10 

.15 

1 

1 

17 

.20 

.01 

.01 

.03 

.10 

.10 

.15 

1 

1 

17 

.50 

.01 

.01 

.03 

.10 

.10 

.15 

3 

3 

19 

1.00 

.00 

.01 

.09 

.09 

.11 

.20 

6 

10 

62 


Table 6. Computational Results for LN3. Solving the network interdiction 
problem within 50% of optimal on network LN3, with 10 units of interdiction 
resoiirce, enumerates 19 near-optimal solutions in .15 seconds. 


The ranges for |A^, |A| and /? are too small to result in appreciable run times for the 
dual-finding and Dijkstra sub-algorithms on these problems. (Of course, the dual-finding 
sub-algorithm is independent of R.) As expected, the run time for the enumeration sub¬ 
algorithm increases dramatically as the number of near-optimal solutions increases and as e 

increases. However, that sub-algorithm is still quite efficient, even for large values of e. 
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As a whole, our interdiction algorithm has proven reasonably efficient, running in 
less than a minute for each of our exartples. It appears that the algorithm would prove an 
efficient tool for somewhat larger networks, and larger values of R, especially on a faster 
conputer. The value of £ will probably be the limiting fector, and this suggests a strategy 
controlling the run time and sheer volume of alternate interdictions by successively 
increasing 8, cautiously, for larger problems. 
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VL CONCLUSIONS AND RECOMMENDATIONS 


This thesis has presented specialized network interdiction algorithms to solve the 
network interdiction problem optimally on planar networks, and also enumerate near- 
optimal solutions. The basic problem is to minimize the maximum flow in a capacitated 
network by destroying arcs using limited resources. We propose two distinct algorithms 
for the problem, the first algorithm being applicable to s-t planar networks (5 and t must lie 
on the perimeter of the network), and the second to general planar networks. 

The first algorithm takes as its input an adversary’s primal network embedded in 
the plane, constructs the dual network, and then utilizes a variant of Dijkstra’s shortest 
path algorithm on an expanded version of the dual to solve the interdiction problem 
optimally. The dual is expanded to deal with the consumption of varying, integral 
amounts of interdiction resource. Using the expanded dual network, the algorithm then 
enumerates all solutions within 100e% of optimal using a modified path enumeration 
algorithm. The algorithm has been implemented and computational results show that the 
methodology is quite efficient. A significant drawback to the algorithm is that it requires 
the network to be s-t planar. We therefore propose an extension of the algorithm to 
handle general planar networks. 

The second, more general algorithm is a fairly straightforward modification of the 
first except that: (a) The algorithm effectively searches for cycles in the dual network 
rather than paths, (b) it must keep track of the number of times (odd or even) that a cycle 
has crossed a specified s-t path, and (c) the shortest path and enumeration sub-algorithms 
must be called once for each arc on that s-t path. Though we do not implement it, the 
algorithm should stiU prove reasonably efficient as it should require (roughly) no more 
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than 2m times more work than the &st algorithm, where m is the minimum number of arcs 
among s-t paths. A suggested modification may further reduce the factor of 2m. 

Further development of the algorithms we have proposed is warranted. The 
second algorithm for general planar networks should be inplemented and tested. An 
efficient planarity testing and embedding algorithm (Hopcroft and Taijan, 1975) should be 
added as a “front end” to these algorithms: Some networks may be planar and susceptible 
to our techniques, but no planar embedding is customarily provided. An automatic 
mechanism for culling duplicate candidate solutions should be appended to our code, also. 
We believe that a sorting procedure employing a hash function would be quite efficient 
and could be applied as candidate solutions are generated (and duplicate solutions never 
presented to the user). Modest modifications to the enumeration algorithm could 
eliminate some of the duplicates, too. 

When a network is non-planar, our methodology is not directly applicable. 
However, if deleting only a few nodes would yield a planar graph, we conjectirre that there 
is some way to extend our methods to fit this situation. The problem of finding the 
minimum number of nodes or arcs whose deletion yields a planar graph is NP-conplete 
(e.g., Garey and Johnson, 1979, p.l95). However, efficient methods exist for finding 
maximal planar subgraphs that are often close to being maximum (Cai, et al., 1993). So, 
identification of nearly planar networks should not be an insurmountable stumbling block 
to extend our methods to such networks. 
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